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A complete analysis of isospin breaking in K amplitudes, including both strong (m„ 7^ m^) and electro- 

magnetic corrections at next-to- leading order in chiral perturbation theory, has been achieved recently pP|. We 
discuss the implication of these effects together with the previously known chiral loop corrections |3I4| . on the 
direct CP- violating ratio e'/e. 



1. INTRODUCTION 

The CP-violating ratio e'/e constitutes a fun- 
damental test for our understanding of flavour- 
changing phenomena within the Standard Model 
framework. The experimental status has been 
clarified by the recent KTEV ,5,, Re (e'/e) = 
(20.7 ± 2.8) • 10"^ and NA48 0, Re (e'/e) = 
(14.7 ± 2.2) • 10""*^, measurements. The present 
world average 



Re (e'/e) = (16.7 ± 1.6) • 10" 



(1) 



provides clear evidence for a non-zero value and, 
therefore, the existence of direct CP violation. 

The CP violating signal is generated through 
the interference of two possible tttt decay 

amplitudes with different weak and strong phases. 
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The isospin amplitudes Aq^2 are defined through 
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In the limit of CP conservation, ^0,^2, and 
are real and positive. In the isospin limit, 
Ai = A2 , X2 = xt i'^ the Standard Model and 
the phases Xi coincide with the corresponding tttt 
phase shifts at i?cm = Af/f. 
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Owing to the well-known "A/ =1/2 rule" , e'/e 
is suppressed by the ratio oj — Rev42/Re^o ~ 
1/22. The strong S-wave rescattering of the two 
final pions generates a large phase-shift difference 
between the two isospin amplitudes, making the 
phases of e' and e nearly equal. Thus, 

TT 
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The large tttt phase-shift difference clearly indi- 
cates that unitarity corrections (final state inter- 
actions) play a crucial role in e'/e |3I4| . More- 
over, this observable is very sensitive to isospin 
breaking effects, because the large ratio l/ui am- 
plifies any potential contribution to A2 from small 
isospin-breaking corrections induced by Aq. 

The CP-conserving amplitudes HeAj , their ra- 
tio u and e are usually set to their experimentally 
determined values. A theoretical calculation is 
then only needed for the quantities ImAj. 

2. THEORETICAL FRAMEWORK 

To obtain the Standard Model prediction for 
e'/e, one starts at the electroweak scale where 
the flavour-changing process, in terms of quarks 
and gauge bosons, can be analyzed in a rather 
straightforward way. Owing to the presence of 
very different mass scales (M^r < Mk <C Mw), 
the gluonic corrections are amplified by large log- 
arithms. The short-distance logarithmic correc- 
tions can be summed up using the Operator Prod- 
uct Expansion (OPE) and the renormalization 
group, all the way down to scales /i < rric. One 
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Figure 1. Evolution from Mw to A'% 0. 



gets in this way an effective AS* = 1 Lagrangian, 
defined in the three-flavour theory |l(Jlllj , 



10 

1=1 



(5) 



which is a sum of local four-fermion operators 
Qi, constructed with the light degrees of freedom, 
modulated by Wilson coefficients Ci (fi) which are 
functions of the heavy masses (M > /i) and CKM 
parameters: 



ait,) = - v,,v;jv^,v:. 



(6) 



Only the yi components are needed to determine 
the CP-violating decay amplitudes. The overall 
renormalization scale fi separates the short- and 
long-distance contributions, which are contained 
in Ci{^) and Qi, respectively. The physical am- 
plitudes are of course independent of /i. 

The Wilson coefiicients are known at the next- 
to-leading logarithmic order |12I13| . This in- 
cludes all corrections of 0«r) and O^+^r), 
where t = ln(Afi/M2) refers to the logarithm 
of any ratio of heavy mass scales Mi,M2 > /i. 
Moreover, the full rrit/Mw dependence (at low- 
est order in a,) is taken into account. 



In order to predict physical amplitudes, one is 
still confronted with the calculation of hadronic 
matrix elements of the four-quark operators. 
This is a very difficult problem, which so far re- 
mains unsolved. Those matrix elements are usu- 
ally parameterized in terms of the so-called bag 
parameters Bi, which measure them in units of 
their vacuum insertion approximation values. 

To a very good approximation, the Standard 
Model prediction for e'/e can be written (up to 
global factors) as ;14.15) 



(1/2) 



(l-f7/B)-0.4B, 



(3/2) 



(7) 



Thus, only two operators are numerically rel- 
evant: the QCD penguin operator Qe governs 
ImAo (A/ = 1/2), while ImAa (A/ = 3/2) 
is dominated by the electroweak penguin oper- 
ator Qg. The parameter fliB takes into account 
isospin breaking corrections, which get enhanced 
by the factor l/u. The value Qjb = 0.25 |16I17| 
was adopted in many calculations |14I18I19] . To- 
gether with Bi ^ I, this produces a large numer- 
ical cancellation in eq. Q leading to unphysical 
low values of s'/e around 7 x 10"'' |14I18| . The 
true Standard Model prediction is then very sen- 
sitive to the precise values of these parameters. 

3. CHIRAL PERTURBATION THEORY 

Below the resonance region one can use global 
symmetry considerations to define another effec- 
tive field theory (EFT) in terms of the QCD Gold- 
stone bosons (tt, K, r]). The chiral perturba- 
tion theory (^PT) formulation of the Standard 
Model j2()l21l22j describes the pseudoscalar-octet 
dynamics, through a perturbative expansion in 
powers of momenta and quark masses over the 
chiral symmetry breaking scale ~ 1 GeV. 

Chiral symmetry fixes the allowed operators. 
At lowest order in the chiral expansion, the most 
general effective bosonic Lagrangian with the 
same SU(S)l ® SU{3)ii transformation proper- 
ties as the short-distance Lagrangian |5j contains 
three terms, transforming as (8l,1_r), (27^,1^;) 
and (8l,8/j), respectively. Their corresponding 
chiral couplings are denoted by gs, 327 and gsw- 

The tree-level K — > tttt amplitudes generated 
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by the lowest-order xPT Lagrangian, 

Gf . 
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98 + 1 927^ [Ml - Ml) - ^fygsw 
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(8) 



do not contain any strong phases. From the mea- 
sured decay rates one gets |23 Iffsl ~ 5.1 and 
I (727 1 ~ 0.29. The gEW term is the low-energy 
realization of the electroweak penguin operator. 

The only remaining problem is the calculation 
of the chiral couplings from the effective short- 
distance Lagrangian jSJ, which requires to per- 
form the matching between the two EFTs. This 
can be easily done in the large-A^c limit of QCD 
|24I25 |. because in this limit the four-quark op- 
erators factorize into currents which have well- 
known chiral realizations: 

98 = - L-2 - - 6i + 64 - I6L5 — -p — j Ce, 



(9) 



Together with eqs. (|SJl, these results are equiv- 
alent to the standard large-iVc evaluations of the 
Bi factors. In particular, for e'/e where only the 
imaginary part of the gi couplings matter eqs. © 



amount to Bg^^^^ 



Bq^^^^ = 1. Therefore, up to 



minor variations on some input parameters, the 
corresponding e'/e prediction, obtained at lowest 
order in both the 1/Nc and xPT expansions, re- 
produces the published results of the Munich 
and Rome ^H] groups. 

The large-iVc limit is only applied to the 
matching between the 3-flavour quark theory and 
xPT, as indicated in Figure ^ The evolution 
from the electroweak scale down to < rric has 
to be done without any unnecessary expansion 
in powers of I/Nq] otherwise, one would miss 
large corrections of the form In {M/m), with 
M ^ m two widely separated scales [^Hl- Thus, 



the Wilson coefficients contain the full fi depen- 
dence. 

The large-A^c factorization of the four-quark 
operators Qi {i 7^ 6, 8) does not provide any scale 
dependence, because their anomalous dimensions 
vanish when Nc — s- 00 26 . To achieve a reliable 
expansion in powers of 1/Nc, one needs to go to 
the next order where this physics is captured |27| . 
This is the reason why the study of the A/ =1/2 
rule has proved to be so difficult. Fortunately, 
these operators are numerically irrelevant in the 
e'/e prediction. 

The only anomalous dimensions which survive 
when Nc — *■ 00 are precisely the ones correspond- 
ing to Qe and Qg jl7l25) . These operators fac- 
torize into colour-singlet scalar and pseudoscalar 
currents, which are /i dependent. This generates 
the factors 



(<?<?) (a*) = - 



f M2 



{mu + md){^i) (ms + md){^J.) 



in eqs. ©, which exactly cancel the fi dependence 
of C6,8(m) at large Nc |17li^6l^7li^?^| . It remains 
of course a dependence at next-to-leading order. 
Thus, while there are large 1/Nc corrections to 
Re{gi), the large- A^c limit can be expected to 
give a good estimate of lTn{gi) [77j . 

4. CHIRAL CORRECTIONS 

The strong phases xi originate in the final 
rescattering of the two pions and, therefore, are 
generated by chiral loops which are of higher or- 
der in both the momentum and 1/Nc expansions. 
Analyticity and unitarity require the presence of 
a corresponding dispersive effect in the moduli of 
the isospin amplitudes. Since the S-wave strong 
phases are quite large, specially in the isospin- 
zero case, one should expect large higher-order 
unitarity corrections. 

The one-loop analyses of -ftT ^ 27r |ll3l4l29j 
show in fact that pion loop diagrams provide an 
important enhancement of the Aq amplitude, im- 
plying a sizeable reduction (~ 30%) of the fitted 
1 58 1 value. This chiral loop correction destroys the 
accidental numerical cancellation in eq. Q, gen- 
erating a sizeable enhancement of the e'/e predic- 
tion 0]. The large one-loop correction to Aq is 
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associated with large infrared logarithms involv- 
ing the light pion mass. 

A complete one-loop calculation, including 
electromagnetic and isospin violation corrections, 
has been achieved recently Q]. It involves the 
0(/) strong frn^ and ©(e^p^) electro- 
magnetic chiral lagrangians. together with the 
non-lep tonic 0{GfP^) :33,34 and 0(Gi?eV) 
|35I36| electroweak lagrangians. 

4.1. 0(p4) xPT 

It is convenient to decompose the isospin am- 
plitudes Ai = Aj e^^' in their different SU{3)l <8) 
SU{3)]i components. The 0(p*) correction to a 
given lowest-order amplitude a\' 



A- 



(X) 



-,ix) 



1 + Ai4^^ + AcAf , (10) 



{X) 



contains a one-loop contribution A l Aj which 
is completely fixed by chiral symmetry plus a lo- 
cal contribution generated by the corresponding 
higher-order chiral lagrangian. The most relevant 
loop corrections take the values )1I3| : 



I^lA^^ = (0.27 ± 0.05) -f 0.47i, 
^lA^^^ = (1.02 ± 0.60) + 0.47 
AlX™^ = (0-27 ± 0.05) -f 0.47i, 
^lA^^"^^ = (-0.04 ± 0.05) - 0.2H, 
Ai4™^ 



(11) 



(-0.50 ± 0.20) - 0.2H. 



The dispersive components depend on the chiral 
renormalization scale , which has been fixed at 
= 0.77 GeV. The quoted uncertainties reflect 
the changes under a variation of between 0.6 
and 1 GeV plus a small contribution from vary- 
ing the short-distance scale /i between 0.77 and 
1.3 GeV. Notice that the most relevant correction 
Ai has a very small uncertainty because it 
is dominated by the non-polynomial part, which 
is associated with the large isoscalar absorptive 
contribution and does not depend on v^. 

(X) 

The local contributions Ac A\ have been 
computed at leading order in the 1/Nc expan- 
sion. At this order, there is matching ambiguity 
because we do not know at which value of the 
chiral scale the estimates apply. This is taken 
into account by the i^^ uncertainty incorporated 



in AlA\^\ Whenever the absorptive loop cor- 

(X) 

rection is large, the final prediction for A} is 
quite insensitive to the values of the low-energy 
chiral couplings adopted in Ac A\^\ 

The CP-conserving parts of the low-energy cou- 
plings Re{gs) and Re(g27) have been fitted to 
the data, together with the phase-shift difference 
Xo - X2- 

4.2. Isospin breaking in e' /e 

There are two sources of isospin breaking: the 
light quark mass difference m,j — and elec- 
tromagnetic corrections. Accounting for isospin 
violation via the general parametrization , the 
ratio u! differs from = Re^J /ReAo by a pure 
AI = 5/2 effect: 



= ^+ (l + /5/2) ; /5/2 = 



ReA2 
ReAt 



1. (12) 



Since is directly related to branching ratios, 
it proves useful to keep tj+ in the normalization 
of e' 37 . The formula for e' takes then the form: 
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where 



Ao 



Im^o Re^, 



(0) 



ReA, 



(0) 



IB 



ReA. 



(0) 



Im^ 



(0) 



(13) 



(14) 



(15) 



and the superscript (0) denotes the isospin limit. 
The quantity 



rioff = f^IB — Ao — /5/2 



(16) 



includes all effects to leading order in isospin 
breaking and it generalizes the more traditional 
parameter fliB- We have adopted the usual (but 
scheme dependent) separation of the electroweak 
penguin contribution to Im^2, ImA™'', from the 
effects of the other four-quark operators. 
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Table 1 

Isospin violating corrections for e' in units of 10^^ 
[2] . LO and NLO denote leading and leading plus 
next-to-leading orders in xPT. 
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^0 




LO 
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LO 


NLO 
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16 ±5 


18 ± 7 


22.7 ±7.6 


Ao 
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-0.4 ±0.1 


9±3 


8.3 ±3.6 


h/2 











8.3 ±2.4 




12 


16 ±5 


9±6 


6.0 ±7.7 



Although r^iB is enhanced by the ratio l/u}^^\ 
the numerical analysis shows all three terms in 
H16|l to be relevant when both strong and elec- 
tromagnetic isospin violation are included. The 
different corrections are shown in Tabled where 
the first two columns refer to strong isospin viola- 
tion only (m,j ^ rrid) and the last two contain the 
complete results including electromagnetic cor- 
rections. Taking a = 0, the isospin breaking is 
completely dominated by the tt'^-jj mixing con- 
tribution OJ'b'' = 0.16 ± 0.03 [HHI. Electro- 
magnetic effects give sizeable contributions to all 
three terms, generating a destructive interference 
and a smaller final value |2I 

f^cff = (6.0 ± 7.7) • 10^2 (17) 
for the overall measure of isospin violation in e'. 

5. DISCUSSION 

The infrared effect of chiral loops generates an 
important enhancement of the isoscalar K — > tttt 
amplitude. This effect gets amplified in the pre- 
diction of e'/e, because at lowest order (in both 
1/Nc and the chiral expansion) there is an acci- 
dental numerical cancellation between the / = 
and 1 = 2 contributions. Since the chiral loop 
corrections destroy this cancellation, the final re- 
sult for e'/e is dominated by the isoscalar ampli- 
tude. The small value obtained for [2 rein- 
forces the dominance of the gluonic penguin oper- 
ator Qq. Taking this into account, the Standard 



Model prediction for e'/e [3] turns out to be 

Re (e'/e) = (l.B ± 0.2 ± 0.5) • 10-^ (18) 

in excellent agreement with the experimental 
measurement The first error has been es- 

timated by varying the renormalization scale /z 
between Mp and TOc- The uncertainty induced 
by rUs ISHl, which has been taken in the range |3] 
{ms + md)(l GeV) = 156 ± 25 MeV, is indicated 
by the second error. 

The most critical step is the matching between 
the short and long-distance descriptions, which 
has been done at leading order in 1/Nc- Since 
all next-to-leading ultraviolet and infrared log- 
arithms have been taken into account, our ed- 
ucated guess for the theoretical uncertainty as- 
sociated with subleading contributions is ~ 30% 
(third error). While a better determination of nis 
can be expected soon, the control of these non- 
logarithmic corrections at the next-to-leading or- 
der in 1 /Nc remains a challenge for future inves- 
tigations 00] ■ 
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